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Abstract 

The photovoltaic array sun tracking control system of 
* Space Station Freedom is described in this paper. A synthesis 
procedure for determining optimized values of the design 
variables of the control system is developed using a 
constrained optimization technique. The synthesis is 
performed to provide a given level of stability margin, to 
achieve the most responsive tracking performance, and to 
meet other design requirements. Performance of the baseline 
design, which is synthesized using predicted structural 
characteristics, is discussed and the sensitivity of the stability 
margin is examined for variations of the frequencies, mode 
shapes and damping ratios of dominant structural modes. The 
design provides enough robustness to tolerate a sizeable error 
in the predicted modal parameters. The paper concludes with 
an investigation on the sensitivity of performance indicators as 
the modal parameters of the dominant modes vary. The 
design variables are re-synthesized for varying modal 
parameters in order to achieve the most responsive tracking 
performance while satisfying the design requirements. This 
procedure of re-optimizing design parameters would be useful 
in improving the control system performance if accurate modal 
data are provided through an on-orbit modal identification 
experiment. 


structure to the inboard structure. The attitude of the inboard 
structure is controlled by Control Moment Gyros (CMG’s) and 
Reaction Control System (RCS) jets. The Alpha Joints are 
used to orient the array surface normal vectors along the solar 
vector so that maximum solar energy falls on the arrays during 
the daylight portion of each orbit. The Alpha Joint control is 
designed to be a basic position tracking system with minor- 
loop velocity feedback to stabilize and provide damping to the 
rigid body tracking motion. A proportional-integral (P-I) 
compensation is added in both the velocity and position loops 
to minimize steady-state tracking error. * 


Port 

side 



Starboard 
side 


EPS Radiator 

X (Flight direction) 
PV Array z (Nadjr) 

Fig. 1 SSF Assembly Complete configuration 


Introduction 

To obtain electric power, Space Station Freedom (SSF) 
shown in Fig. 1 depends on photovoltaic (PV) solar arrays 
which track the sun during orbital daylight. The arrays are 
attached to deployable masts which are in turn attached 
through a rotary joint, called a Solar Beta Rotary Joint (Beta 
Joint), to the outboard portion of the transverse booms. The 
Beta Joint permits rotation of the arrays to compensate for the 
variation of the orbit plane with respect to the ecliptic plane. 
Rotary joints, called Solar Alpha Rotary Joints (SARJ or Alpha 
Joints), regulate the relative rotational position of the outboard 
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The allowable rigid-body control bandwidth of the Alpha 
Joint controller encompasses the resonant structural 
frequencies of the outboard boom and photovoltaic system so 
that the possibility for adverse interaction between the control 
system and the dynamic response of the structure exists. To 
reduce the possible detrimental effect of control/structure 
interaction, a low-pass filter is added to the velocity loop to 
attenuate the structural response signal. The proper 
placement of the comer frequency of the filter and selection of 
values for gain setting of the P-I compensation in the velocity 
and position loops are required to provide optimum 
performance. The proper selection of these design values 
depends on the accuracy of the predictions of the structural 
frequencies and modal response at sensor locations. 

The space station is too large and flexible to support its 
own weight on earth. Hence, the structural dynamic 


characteristics will have to be estimated from analytical models 
and component modal tests rather than from modal tests of the 
assembled structure. Further, the station is designed to support 
reconfigurations for more power, payload installations and other 
activities which can change the dynamic characteristics of the 
station. Because of the considerable uncertainty involved in 
predicting the dynamic characteristics of the station initially and 
over its lifetime, the sun tracking control system must be 
designed wilh a high degree of stability robustness. This paper 
presents results from an investigation of the sensitivity of the 
control system stability margins to variations in modal 
parameters of dominant structural vibration modes. Wilh this 
sensitivity established, the control system can be designed with a 
reasonable degree of robustness to assure stable tracking for a 
given range of variation in structural parameters which might 
occur due to configuration changes and errors in analytical 
estimation. 

The paper describes the procedures to attenuate the 
possible control/struclurc interaction and to examine controller 
sensitivity to variations in structural modal parameters. First, the 
significant components of the space station related to the Alpha 
Joint control are described and the sun tracking control system is 
described. A baseline design is then determined using 
constrained optimization techniques to meet design requirements, 
to provide a given level of stability margin, and to obtain the most 
responsive tracking possible consistent with the assumed 
structural characteristics. Performance of the baseline design is 
discussed and the sensitivity of the stability margin is examined 
for variations of the natural frequencies, mode shape amplitudes 
and modal damping ratios of the dominant structural modes. The 
paper concludes with an investigation on the sensitivity of 
performance indicators as the modal parameters of the dominant 
modes vary. The design variables are re-synthesized for varying 
modal parameters in order to achieve the most responsive 
tracking performance while satisfying the design requirements. 

Description of space station and 

sun tracking function 

Description of space s tation structure 

The space station structure can be broadly divided into an 
inboard core structure and an outboard articulating structure. As 
shown in Fig. 1, the inboard core structure is comprised of a 
module cluster, center truss, thermal control system (TCS) 
radiators, and various user payloads. The port and starboard 
truss, PV arrays and electrical power system (EPS) radiators 
constitute the outboard articulating structure commonly referred 
to as the outboard structure. The attitude of the inboard core is 
maintained close to a local-vertical-local-horizontal (LVLH) 
orientation using active control devices such as the CMGs and/or 
RCS thrusters. The LVLH X-axis is parallel to the flight 
direction, the LVLH Z-axis directed to tire nadir, and the LVLH 
Y-axis is orthogonal to the orbit plane. The power required for 
the space station is provided by the PV arrays as they track the 
sun. Due to the motion of the space station along its orbit, core 
attitude fluctuations with respect to the LVLH orientation, and 
ihe variations of the orbit geometry wilh respect to the sun, the 


orientation of the PV arrays have to be constantly adjusted with 
respect to the inboard core to track the sun. This function is 
performed by the Alpha Joints and the Beta Joints. The Alpha 
Joints provide a relative rotational motion between the inboard 
core structure and the outboard articulating structure. The Beta 
Joints perform the PV array orientation adjustment with respect 
to the articulating outboard port and starboard truss. 

Assuming that the station is maintained at an LVLH 
attitude, the Alpha Joint rotation rate would be the orbital rate, 
completing a revolution every 90 minutes. The Beta Joint 
rotation is extremely slow over an orbit and follows closely the 
yearly variation of the orbit plane. Therefore, the Beta Joint 
drive and control is not addressed in this paper. 

Physical -dcscriptiQn oL Alpha Joint 

Each Alpha Joint consists of dual motors, dual resolvers, a 
motor controller, drive pinions, a bull gear and trundle bearings 
as depicted schematically in Fig, 2,^ The motor provides the 
torque required to rotate the outboard structure. The amount of 
control torque provided by the motor is determined by the motor 
controller based on the measurements obtained by the joint 
resolver and the desired rotation data from the Velocity Vector 
Generator (VVG) located on the inboard structure. The motor 
drive pinion to bull gear ratio has been carefully selected to 
minimize mechanical parts count and hence maximize 
reliability.** The bull gear is rigidly attached to the outboard 
structure through a shear plate. The large bull gear (about 10 
feet in diameter) is equipped with trundle bearings in order to 



Fig. 2 Schematic of Alpha Joint drive train and control system 
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accommodate large temperature gradients. The trundle bearings 
arc the main source of friction. A set of high power roll rings 
(not shown in Fig. 2) carries power across the joint as the joint 
rotates.** 

Description of Alpha joint control system 

The SARJ motor controller generates the required motor 
torque based on the difference between the desired and measured 
relative joint angular position and angular velocity at points A 
and B as shown in Fig. 2. Point A is fixed to the inboard structure 
while point B is located on the outboard structure. The desired 
relative joint velocity command is determined by the VVG 
located on the inboard structure and is an input to the SARJ motor 
control system. As the SARJ position leads or lags the desired 
position (also provided by the WG), the input to the velocity- 
loop is increased or decreased to compensate for the position 
error. Hence the reference input to the control system includes 
the desired relative angular velocity (co B - co A ) and the desired 
relative angular position (0 B - 0 A ). 

The detailed control system used in this paper is based on 
the design obtained from the SSF Preliminary Design Review 
document. 4 A block diagram of the control system is shown in 
Fig. 3. The control system consists of an inner velocity servo 
loop which tracks the desired joint velocity and an outer position 
servo loop which increases or decreases the velocity command 
depending upon the position error. The velocity command is 
converted to a voltage command and the maximum allowable 
SARJ velocity is maintained by a voltage limiter. 

The inner servo loop includes a fourth order Butterworth 
bending filter^ with two additional zeros to roll-off the high 
frequency structural modes. The two additional zeros are 
included to reduce the loss of rigid-body phase margin due to 
phase shift. The transfer function of the bending filter is 


<*■>■< 4 <0 

f s 4 + 2.6o c s 3 + 3.4cd 2 s 2 + 2.6to 3 s + co 4 

where the comer frequency (co c ) and the frequency of the zeros 
(co 0 ) arc design parameters. 

While the inner velocity feedback loop increases system 
damping, it also increases steady-state tracking error. Integral 
control helps to reduce steady state errors. Therefore, a 
proportional-integral (P-I) controller is used for the inner 
velocity loop and a similar P-I controller is used for the outer 
position loop. The position loop has a double integrator (one in 
the outer loop and the other in the inner loop) to track a ramp 
signal with zero steady state error. The P-I controllers are also 
provided with integration limits to prevent the system from being 
overdriven. 4 This is required due to the acceleration limits 
imposed on the SARJ. Each P-I controller has two gain settings: 
kp and kj are the proportional and integral gains for the velocity 
loop and k' p and k\ are the corresponding gains for the position 

loop. These are also included as design variables during the 
synthesis of the control system design. 

The power amplifier shown in Fig. 3 is accompanied by a 
current limiter. A 66 in-lbf/amp motor torque constant and a 
gear ratio of 325 is assumed. 4 The output torque from the gear is 
subject to the large static and dynamic friction of the trundle 
bearings. The output torque must exceed the static friction of 
3580 in-Ibf to initiate motion of the SARJ. Once the motion is 
initiated, a net torque, which exceeds the dynamic friction of 
2870 in-lbf, is applied to the structure at point B (actuator point) 
as shown in Fig. 2. 

A finite-element model of the space station structure is 
used to compute the natural frequencies and mode shapes. A 
flexible body state space model of the structure is formed from 
the modal data. This is combined with the rigid body model of 
the outboard structure to characterize the dynamics of the space 
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slaiion structure. The rigid body inertia of of the outboard 
structure (cither port or starboard side) about the Y-axis (Iyy) * s 

1.75xt0 7 lbf-in-sec 2 . Details of the finite -clement model and 
modal description arc given in the next section. The net torque 
(T n ) applied to the structure at actuator point B (Fig. 2) causes 

motion of the flexible structure. The net torque applied on the 
outboard structure causes a reaction torque on the core structure. 
This reaction torque has to be compensated by RCS jet and/or 
CMG torques to maintain the attitude of the core structure. This 
paper assumes that the reaction torque is compensated ideally and 
the rigid body attitude of the core structure remains stationary. 
The joint position and velocity containing both rigid body and 
elastic components arc measured and fed back to close the control 
system loops. 



Fig. 5 Undamped natural frequencies of SSF below 5 Hz 


Dynamic characteristics of Space Station Freedom 

A finite element model of SSF was used to investigate the 
influence of clastic response to the Alpha Joint control system 
performance. To provide a rotational dcgrec-of- freedom (DOF) 
about the Y-axis of the outboard structure, two coincident grids 
were placed at the center node of the Alpha Joint as indicated in 
Fig. 4. These grids are rigidly connected in the other Five DOF. 
The undamped natural frequencies of the finite -element model 
arc shown in Fig. 5. For the model used, there are 240 modes 
below 5 Hz including the eight rigid body modes. 



Fig. 4 Finite element model of Space Station Freedom 


Since the outboard structure is subject to a continuous 
rotational motion while the SARJ is actively controlled, the 
influence of the outboard structure orientation on the SSF 
dynamic characteristics was investigated. The outboard structure 
was rotated 90 degrees from the minimum drag PV orientation 
shown in Fig. 4 and the natural frequencies and mode shapes were 
calculated. The resulting maximum drag PV orientation SSF 
exhibits a natural frequency distribution almost identical to that 
of the minimum drag orientation. Also, the transfer function 
characteristics between the two points where the SARJ controller 
is applied (points A and B in Fig. 2) were almost identical. Thus 
the dynamics of SSF with respect to the SARJ control is 
considered to be lime invariant for this study. Moreover, due to 
the symmetry of the space station structure, the dynamic 
characteristics of the SARJ at the port and starboard sides are 
nearly identical. Thus, as a representative plant model, the port 
SARJ with a minimum drag PV array orientation is employed in 
this paper. 


The equation of motion governing the flexible response at 
the port SARJ is represented by 

x = Ax + bu n \ 

y = Cx K) 


where u is the control torque and 



where q is the modal displacement vector; 0 A and 0 B are the 
angular displacements at the points A and B, respectively; co A 
and g> b are the angular velocities at the points A and B, 

respectively; Q = diag{o).} and Z = diag{£.) in which ©. and Q 

are the natural frequency and modal damping ratio of the ith 
mode, respectively; <|> A and <f> B are the row vectors of the unity 

mass normalized mode shape matrix corresponding to the Y- 
rotational DOF at the points A and B, respectively. Modal 
damping ratio of 0.1% for all the flexible modes is assumed as a 
baseline value for the design and simulation of the SARJ control 
system. 

Synthesis procedure of control system 
design variables 

CtiniroLsYstenr requirements 

The SARJ control system discussed earlier has six design 
variables whose values can be adjusted to optimize control system 
performance and to satisfy prescribed requirements. The control 
system requirements can be classified as frequency-domain and 
time-domain requirements. The frequency domain requirements 
in both inner and outer loops are as follows: (1) Rigid-body 
open-loop gain and phase margins (GM and PM) must be greater 
than or equal to 6 dB and 45 degrees, respectively, to assure a 
stable rigid-body motion; (2) Apparent gain margin (AGM), 
defined as the minimum distance of the open-loop gain from the 
zero dB line in the frequency range encompassing the structural 
resonance frequencies, should be at least 20 dB to guarantee 
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sufficient stability robustness to uncertainties in the predicted 
modal parameters; and (3) The closed-loop poles associated wilh 
rigid body and controller should have a minimum damping ratio 
of 0.5. By constraining rigid body and controller closed-loop 
poles to a prescribed sector in the complex plane, this frequency- 
domain constraint assures low overshoot during the transient 
response. The time-domain requirements 4 include small steady- 
state pointing errors, small peak-to-pcak variation of steady-state 
pointing errors over long lime intervals and low jitter. The jitter 
is defined as the peak-to-peak variation of the position error in 
one second. 

The Preliminary Design Review document 4 requires that 
the inner and outer closed-loop bandwidths (BW V and BWp) be 

between 0.01 to 1 Hz. This requirement is treated in this paper as 
part of a performance index which is to be maximized. 
Bandwidth is a measure of the responsiveness of a control system 
(closely related to the rise time) and also represents a disturbance 
rejection threshold. Core attitude fluctuations and the maneuvers 
for feathering and debris collision avoidance may demand large 
attitude rate changes and thus motivate higher responsiveness of 
the control system. The other component of the performance 
index to be maximized is the magnitude of the real part of the 
dominant rigid-body and controller closed-loop pole (a). The 
dominant pole is defined here as the rigid-body and controller 
closed-loop pole (for both inner and outer loops) closest to the 
imaginary axis. This is equivalent to minimizing a settling time in 
a time-domain analysis. Table 1 summarizes the design objective, 
the control system frequency-domain and time-domain 
requirements, and the design variables. 


Table 1 Summary of Alpha Joint controller design objective, 
variables, and constraints 


Design 

objective 

Maximize position and velocity loop bandwidth 
and minimize settling time 

Design 

variables 

Controller gains (kp, k k p and k '; ) 

Compensation filter break frequencies (to c and( 0 o ) 

Constraints 

Rigid-body gain margins ^6 dB 
Rigid-body phase margins >45 deg 
Apparent gain margins in structural resonant 
frequency range £ 20 dB 
Minimum rigid body and controller 
damping ratio ^0.5 
Bandwidth: between 0.01 and 1 Hz 
Steady -state pointing error ^0.58 deg 
Peak-to-peak variation of steady-stale pointing 
error over 30 min: between ±0.5 deg 
Jitter <0.01 deg/sec 


Synthesis procedure 

The plant model includes the rigid body inertia of the 
outboard structure about the Alpha Joint axis and the flexible 
modes of the entire structure. Figure 6 shows a frequency 
response function (FRF) of the rigid plant only compared with 
the FRF of the rigid body wilh all flexible modes of the Finite 
clement model up to a frequency of 5 Hz. The FRF shown is the 
magnitude ratio of the velocity response to net torque applied and 
is plotted using a decibel scale against the log of the frequency. A 
control system designed with only the rigid plant taken into 


account could become unstable if the dominant flexible modes are 
not well attenuated by the control system. The dominant flexible 
modes (or dominant modes) are the most influential modes 
among the flexible modes in determining the apparent gain 
margin. The flexible modes are most likely to be influential if 
their gain is high and their frequency is low. Figure 6 indicates 
that the dominant flexible modes occur at frequencies of 0.485 Hz 
and 0.486 Hz. The corresponding mode shapes are shown in Fig. 
7. These modes correspond to a rigid body rotation of the 
outboard trusses coupled with symmetric and anti-symmetric 
bending of the PV arrays while most parts of the inboard 
structure remain stationary. Other modes which might interfere 
with rigid body controllers are at higher frequencies. Even 
though their gains might be higher than the modes selected as 
dominant, their influence would be attenuated further by a low- 
pass filter used to roll-off the effects of the dominant modes and 
can hence be ignored during the control system design. 



Frequency, Hz 


Fig. 6 Frequency response of rigid body and rigid+fiex body 
plants (velocity response to control torque) 



Fig. 7 Dominant modes for SSF Alpha Joint control 


It is desirable that all the flexible modes be removed from 
the plant during the synthesis of design variables to ease the 
computational load. The dominant flexible modes have a gain of 
approximately 65.4 dB above the rigid body gain as shown in Fig. 
6. Hence the 20 dB apparent gain margin constraint for the plant 
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with the flexible modes is equivalent to a constraint of 85.4 dB 
apparent gain margin for the rigid body plant at the frequency of 
0.49 Hz. For the synthesis of design variables, the plant is 
coasidered as a rigid body with one of the constraints modified as 
described above. This consideration is used only for design 
purposes and not for subsequent frequency response analyses or 
lime response simulations. 

To further simplify the synthesis procedure, the limiters in 
the control system are ignored. The friction block is replaced by 
a transfer function of unity, which is a conservative assumption 
for a robust design. The linearized plant and control system are 
transformed to the frequency domain and the design synthesis is 
performed using a coastrained optimization scheme. 

Now, the synthesis problem can be stated as follows: Find 
the optimum values for the six design variables 

(co c , 0 ) o , kp, k., kjj, and k!) which maximize the performance 
index 

J = p 1 BW v + M 2 BW p +p 3 o (3) 

while satisfying the following constraints: 

Velocity -loop rigid-body gain margin > 6 dB 
Position-loop rigid-body gain margin > 6 dB 
Velocity-loop rigid-body phase margin > 45 degrees 
Position-loop rigid-body phase margin > 45 degrees 
Velocity-loop rigid-body gain 

at the frequency of 0.49 Hz < -85.4 dB 
Minimum rigid-body position and velocity 
closed-loop damping factor > 0.5. 

The scalars p.p and p. 3 are weighting factors selected to give 

equal weights to each element of the performance index. Equal 
weights of unity are used for Pj and p 2 since velocity and 

position closed-loop bandwidths are of the same magnitude. The 
magnitude of a is an order of magnitude less than the bandwidths 
expressed in rad/sec. To give approximately equal importance to 
the settling time, p 3 = 10 is selected. The control torque is not 


incorporated in the synthesis problem either as a performance 
index or as a constraint. It is assumed that the available control 
torque will be sufficient to perform the maneuver required by the 
synthesized design. 

The constrained optimization problem is solved using a 
nonlinear programming method available in MATRIXx 
software.** A brief description of the method is given here. The 
constrained nonlinear optimization problem is initially 
approximated by a linearly constrained optimization problem 
with an augmented Lagrangian objective function. Then, 
sequential quadratic programming is implemented to solve the 
optimization problem with the new objective function 
approximated to the second order. The software uses a recent 
extension^ of the Karmarkar’s interior point algorithm to solve 
the resulting quadratic programming problem. It was noticed 
that there existed more than one local minimum for the 
optimization problem so that selection of the initial guess values 
for the six design variables was important in obtaining the "best" 
local minimum. However, one cannot be guaranteed that this 
"best" local minimum obtained is the global minimum. 

The resulting optimized values of the design variables with 
performance index and constraints are summarized in Table 2. 
The design satisfies all the prescribed design requirements. The 
steady -state time domain requirements listed in Table 1, which 
were not enforced during the optimization, were checked for 
violation using the optimized design variables through time 
response simulation. The closed-loop poles of the velocity and 
position loops were inspected to ensure stability. The following 
section discusses the design results in detail. 

Design results 

The results of the design are summarized in Table 2. The 
bandwidths of the velocity and position loops are within the range 
specified in the space station program requirements. Design 
constraints are all met near or at the boundary of the constraints. 
The compensated Bode plot of the velocity loop on the port side 


Table 2 Results of Alpha Joint controller design 



Description 

Requirements 

Optimized results 


Velocity-loop bandwidth (BWv) 

0.01 Hz < BWv ^1 Hz 

0.053 Hz 

Design 

Position-loop bandwidth (BWp) 

0.01 Hz <BWp < 1 Hz 

0.027 Hz 

objectives 

Distance of dominant pole to imaginary axis 

N/A 

0.026 


CD c 

N/A 

0.54 rad/sec 


C0 o 

N/A 

3.83 rad/sec 

Design 

Xp 

N/A 

0.73 

variables 

k i 

N/A 

0.016 


k P 

N/A 

1.07 



N/A 

0.023 


Rigid-body velocity-loop gain margin 

£6dB 

7.5 dB 


Rigid-body position-loop gain margin 

>6dB 

7.6 dB 

Design 

Rigid-body velocity-loop phase margin 

^45 deg 

45.0 deg 

constraints 

Rigid-body position-loop phase margin 

>45 deg 

45.0 deg 


Velocity-loop apparent gain margin at 0.49 Hz 

£20 dB 

20.0 dB 


Minimum rigid-body and controller damping ratio 

£0.5 

0.51 
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of the station is shown in Fig. 8. Ail the flexible modes below 5 
Hz were incorporated in the simulation. An apparent gain 
margin of 20 dB and rigid-body gain and phase margins of 7.5 dB 
and 45 deg, respectively, arc obtained as indicated in the figure. 



Frequency, Hz 

Fig. 8 Compensated Bode plot of port velocity open loop 



Fig. 9 Simulated position and velocity tracking response 


shown in Fig. 10. The steady-state jitter requirement is met in 
less than one minute. 



Fig. 10 Simulated position tracking error and jitter 



Fig. 1 1 Simulated motor/gear torque and net torque applied 


The time response of the control system is simulated for a 
step velocity command of 4 deg/min and a ramp position 
command with a slope of 4 deg/min. Figures 9 through 1 1 show 
the results of the simulation using all the flexible modes below 5 
Hz. The position command input and the resulting response arc 
compared in Fig. 9. No motion occurs until the static friction is 
overcome. Then, after a brief initial transient period, the 
tracking is performed accurately. The velocity response to the 
command is also shown in Fig. 9. The higher frequency 
component of the response corresponds to the flexible response 
of the structure at the dominant mode frequency. The flexible 
response is quite small and should cause no structural load 
problems. The position error, i.c., the difference between the 
position command and the actual position response, is shown in 
Fig. 10. The steady state pointing error, which should be less 
than 0.58 deg, is met within 40 sec. Jitter time history is also 


Figure 1 1 illustrates the torque generated by the motor and 
gear train to perform the Alpha Joint pointing and the net torque 
applied to the station after overcoming the friction. No net 
torque is applied to the station until the torque overcomes the 
static friction. Once the motion of the Alpha Joint is initiated, the 
dynamic friction becomes effective and the magnitude of steady- 
slate motor/gear torque is just enough to overcome the 
dynamic friction. As a result, the core structure of the station 
docs not experience a net torque applied from the joint motor and 
gear train after the initial transient period of 100 sec. The 
maximum torque available is approximately 30,000 in-lbf. The 
peak motor torque required is well within the torque limit. 
Therefore, the assumption that the motor is capable of producing 
Lhe level of torque required by the synthesized design is verified. 
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Sensitivity of stability robustness and performance 
to structural parameter variations 

The flexible plant model used in the synthesis of design 
variables is characterized by its natural frequencies, mode shapes, 
and modal damping ratios. The baseline design variables 
presented in the previous section was obtained using "nominal" 
values of the modal parameters. These nominal modal 
parameters are analytical estimates and may be inaccurate. In 
order to address the effect of variation of the modal parameters 
on the control system performance, two types of sensitivity 
analyses are performed in this paper. The first study examines 
the change in the apparent gain margin (stability robustness 
measure) as the modal parameters vary while maintaining the 
baseline values of the design variables. In the second study, the 
values of the design variables are re-synthesized for the flexible 
plant with modified modal parameters while enforcing the same 
design constraints as used for the baseline design. The resulting 
variations of the performance index and its components are 
investigated. For a set of fixed values of the design variables, the 
performance index is independent of the changes in structural 
modal parameters because the performance index is composed of 
the velocity and position loop bandwidlhs and the dominant rigid- 
body closed-loop pole. Thus, performance robustness due to 
structural parameter variation is not discussed in this paper. 

Sensitivity of stability robustness 

While keeping the baseline values of the design variables 
unchanged, the effect of modal parameter changes on the stability 
margin is examined. The gain and phase margins of the position 
and velocity loops and the locations of the rigid-body plant and 
controller poles remain unchanged since the values of the design 
variables and rigid body plant are fixed for this study. Only the 
apparent gain margin is influenced by the variation of the modal 
parameters. The baseline design parameters obtained for the 
nominal flexible plant provide 20 dB velocity-loop apparent gain 
margin at the frequencies of dominant modes. The position-loop 
always has a larger apparent gain margin. Therefore, this study 
examines only the change in the velocity-loop apparent gain 
margin as the modal parameters of the dominant modes vary. 

Figure 12 shows the sensitivity of stability margin to a 
variation in dominant mode frequencies and modal amplitudes 
while keeping the modal damping ratio constant at 0.1%. Lines 
of constant apparent gain margins are shown in the figure. 
Changes in modal frequencies affect only the system A matrix. 
The modal amplitude is varied by positive scalar multiplications 
of the mass normalized mode shapes at the sensor and actuator 
locations simultaneously. These affect the system b vector and C 
matrix. All modal parameter variations are represented in 
percentages. 

As the dominant mode frequencies decrease, the resonant 
peaks of the two dominant modes move closer to the comer 
frequency of the bending filter (co c ). Hence, the amount of roll- 
off at the dominant mode frequencies decreases resulting in a 
smaller apparent gain margin. This is indicated along the 
ordinate of Fig. 12. When the dominant mode frequencies are 
decreased by approximately 38% while keeping the modal 


amplitude constant, the apparent gain margin diminishes to zero 
and if the phase at this frequency is at -180 degrees the system 
becomes unstable. The origin corresponds to the nominal plant. 
It should also be noted that the robustness would be improved if 
the dominant mode frequencies are increased (not shown in Fig. 
12). For a given frequency change, the reduction in robustness is 
seen along the abscissa of Fig. 12 as the modal amplitude 
increases. The increase in modal amplitude corresponds to the 
increase in height of the resonant peaks and hence reduces the 
apparent gain margin. A change of approximately +210% in the 
modal amplitude of the dominant modes is required to nullify the 
apparent gain margin. A reduction in modal amplitude increases 
the apparent gain margin (not shown in Fig. 12). 



Fig. 12 Sensitivity of velocity-loop apparent gain margin to 
variations in frequency and modal amplitude of dominant 
elastic mode (modal damping = 0.1 %) 


Points interior to the axes correspond to combined errors 
in both modal frequencies and modal amplitudes. Regions of 
negative stability margins are also shown in Fig. 12. Space 
station structural parameters predicted analytically may contain 
appreciable errors.** The inability of ground testing for model 
verification and the synthesis error of component modal 
characteristics are the major source of the errors. Thus, in order 
to compensate for the errors, the control system design should be 
robust to large variations in modal parameters. Figure 12 
indicates that, for a nominal design using 20 dB apparent gain 
margin, the control system can tolerate a wide range of modal 
parameter variations. 

Figure 13 depicts the sensitivity of stability robustness to 
a variation in modal damping ratio while keeping the dominant 
mode frequencies and modal amplitudes constant. The nominal 
plant has a conservative low damping of 0.1%. As the modal 
damping ratio increases, the resonant peak of the dominant modes 
decrease (the resonant peak is approximately proportional to 
1/2Q. Thus, the apparent gain margin increases. As the modal 
damping ratio decreases the apparent gain margin decreases. 
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However, this is not discussed due to the already conservative 
nominal value of modal damping ratio chosen. 



In this study, the control system design variables arc re- 
synthesized as the modal parameters of the dominant modes vary. 
The constraints, as listed in Table 1, are enforced and the values 
of the design variables arc re-computed for each change in' the 
value of a modal parameter to provide the most responsive 
tracking performance. The variation of the performance index 
and its three components is examined as the frequencies, modal 
amplitudes, and modal damping ratios change. 



Fig. 13 Sensitivity of velocity-loop apparent gain margin to 
variations in modal damping ratio of dominant modes 



Change in frequency of dominant modes from 
assumed nominal frequency, % 

Fig. 14a Sensitivity of performance indicators to variation in 
frequency of dominant flexible modes 

Figure 14a shows the variation of the performance 
indicators as the frequencies of the dominant modes change. The 
modal amplitudes and damping ratios are held at their nominal 
values. The ordinate is the ratio of the performance index (J) and 
its indicators (BW V , BW p , and o) to their respective nominal 


values, The performance index consistently improves as the 
dominant mode frequencies increase and vice versa. As the 
dominant mode frequencies increase, the passband of the 
Butterworth filter can be extended. This results in an increase in 
the performance index and its components. The design variables 
arc re-synthesized at the interval of 5% change in nominal 
frequency. For instance, the baseline design variables are used as 
an initial guess for the synthesis with 5% frequency change. 

Figure 14b illustrates the variation of the performance 
indicators as the modal amplitude of the dominant modes 
changes. The frequencies and modal damping ratios are fixed at 
their nominal values. The performance index and its components 
consistently decrease with increasing modal amplitude. As the 
resonant peaks of the dominant modes increase, the apparent gain 
margin constraint is violated. In order to meet the constraint, the 
comer frequency of the filter decreases while ensuring that the 
other constrains are not violated. Thus, the performance index 
and its components, in general, decrease. The design variables 
are re-synthesized at intervals of 10% change in nominal modal 
amplitude using the previous optimized design as a starting guess 
for design variables. 



^ -40 -20 0 20 40 


Change in modal amplitude of dominant modes 
from nominal modal amplitude, % 

Fig. 14b Sensitivity of performance indicators to variadon in 
modal amplitude of dominant flexible modes 

The variations in the performance indicators with respect 
to the changes in modal damping rarios are shown in Fig. 14c. As 
the modal damping ratio increases, the resonant peaks of the 
dominant flexible modes reduces and the apparent gain margin 
constraint becomes less critical. The passband of the low-pass 
filter may be extended, increasing the bandwidth until the 
apparent gain margin constraints or the other constraints become 
active. Hence the performance indicators should increase or at 
worst remain constant. Figure 14c indicates that the performance 
increases until the modal damping ratio is raised to 
approximately 0.4%. At this value, the apparent gain margin 
constraint is no longer active and the performance indicators 
remain constant for additional increase in damping. Thus, 
attempts to increase the modal damping ratio would be beneficial 
in increasing control system performance only up to a certain 
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level. The magnitude of modal damping ratio at which this 
performance saturation occurs will depend on the constraints 
imposed. 

The lime domain analysis for each set of optimized values 
of the design variables should be performed lo check whether the 
steady-state pointing and jitter requirements arc satisfied. The 
time domain response will depend on the values of the design 
variables but it will damp out since the closed-loop rigid-body 
pole constraint ensures a minimum damping ratio of 0.5. The 
steady-state tracking errors could be met by adjusting the 
integrator limits in the controllers. 



Fig. 14c Sensitivity of performance indicators to variation in 
modal damping ratio of dominant flexible modes 


Conclusions 

The significant components of the space station related to 
the Alpha Joint control were described and the photovoltaic array 
sun tracking control system was presented. A baseline design was 
determined using constrained optimization techniques to meet 
design requirements, to provide a given level of stability 
robustness, and to obtain the most responsive tracking capable 
consistent with the assumed structural characteristics. 
Performance of the baseline design was discussed and the 
sensitivity of the stability margin was examined for variations in 
the natural frequencies, mode shapes and damping ratios of 
dominant structural modes. 

The dominant clastic modes occurred approximately at 0.5 
Hz and were well below the maximum Alpha Joint controller 
closed-loop bandwidth of one Hz specified in the design 
requirements. These modes exhibited symmetric and anti- 
symmetric photovoltaic array bending coupled with rigid 
rotation of the transverse booms outboard of the Alpha Joint. 
The highest bandwidth achievable for a design which met all 
requirements for phase and gain stability robustness, steady state 
behavior and reasonable settling time was 0.054 Hz for the 
velocity loop, approximately a decade below the dominant mode 
frequencies. The design provided enough robustness to tolerate a 
sizeable error of up to 40 percent in the predicted resonant 
frequency of the dominant modes and to tolerate an error in the 


predicted modal amplitude of up to 200 percent The station will 
not be tested as a complete system before being placed in orbit 
and modal frequencies will have to be predicted using 
information from component tests and unvalidated analytical 
models. Hence, the high levels of stability margin suggested in 
this study for the nominal design seem appropriate. 

The paper concludes with an investigation on the sensitivity 
of performance indicators as the modal parameters of the 
dominant modes vary. The design variables are re-synthesized 
for varying modal parameters in order to achieve the most 
responsive tracking performance while satisfying the design 
requirements. This procedure of re -optimizing design 
parameters would be useful in improving the control system 
performance if accurate modal data are provided through an on- 
orbit modal identification experiment. 
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